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Modelling the time variation of the surface differential rotation in 
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Sequences of Doppler images of the young, rapidly rotating late-type stars AB Dor and LQ Hya show that their equatorial 
angular velocity and the amplitude of their surface differential rotation vary versus time. Such variations can be modelled 
to obtain information on the intensity of the azimuthal magnetic stresses within stellar convection zones. We introduce a 
simple model in the framework of the mean-field theory and discuss briefly the results of its application to those solar-like 
stars. 



© 2007 WILEY-VCH Verliig GmbH & Co. KGaA, Wcinhcim 



43 

6 : 



> 
in 

(N 
(N 



o 



X 



1 Introduction 

Doppler imaging techniques allow us to measure the sur- 
face differential rotation in rapidly rotating late-type stars 
by tracking the longitudinal motion of starspots located at 
different latitudes (Collier Cameron 2007). Specifically, the 
surface angular velocity at colatitude 9 is assumed to be 
given by a solar-like law: 

Q(9) = Q oq ~dncos 2 9, (1) 

where Q eq is the equatorial angular velocity and d£l is the 
pole-equator angular velocity difference. fl cq and dfl can be 
measured by fitting the shear of starspots along sequences 
of photospheric images covering successive rotations. Al- 
ternatively, fl and d£l can be included as free parameters 
in a code that reproduces the line profile distortions due to 
starspots, thus obtaining their values by a suitable \ 2 mini- 
mization when a sufficiently long time series of line profiles 
is available. 

The long-term monitoring of surface differential rota- 
tion of the two late-type dwarfs AB Doradus and LQ Hy- 
drae shows that their equatorial angular velocity and surface 
shear are functions of the time (see Donati et al. 2003a; Jef- 
fers et al. 2007). It is interesting to note that the variations of 
f2 eq and dQ, are compatible with an internal angular veloc- 
ity uniform on cylindrical surfaces co-axial with the rotation 
axis (Donati et al. 2003a). 



2 Modelling differential rotation variations 

Differential rotation is the result of the re-distribution of an- 
gular momentum in a stellar convection zone under the ac- 
tion of meridional circulation, Reynolds stresses (i.e., the 



correlations of the components of the turbulent velocity in 
a rotating star), and torque by the Lorentz force (quantified 
by the azimuthal Maxwell stresses). The shear associated 
with differential rotation is opposed by the turbulent viscos- 
ity present in the convection zone that tends to make the star 
rotate rigidly. 

In order to model the time variation of the differential 
rotation, let us consider the equation for the angular mo- 
mentum conservation in an inertial reference frame in the 
framework of the mean-field theory (e.g., Riidiger & Holler- 
bach 2004): 



^(pr 2 sin 2 9 V) 



V = 0, 



(2) 



where r is the distance from the centre of the star, t the 
time, p the density, f2(r, 9, t) the angular velocity and the 
angular momentum flux vector given by: 

= (pr 2 sin 2 6T2)it( m ) 

+ r S m9(pu / u' v ) - T -^{BB V + {B'B'J), (3) 
M 

where u m is the meridional circulation, u' the fluctuating 
velocity, B the mean magnetic field, B' the fluctuating 
magnetic field and jl the magnetic permeability; angular 
brackets () indicate the Reynolds mean, defining the mean 
quantities. Note that in Eq. (O the first term specifies the 
angular momentum transport by meridional circulation, the 
second by Reynolds stresses and the third by Maxwell 
stresses, respectively. 

The Reynolds stresses can be written as: 



, , , dui duj , 



(4) 
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where the first term represents the effect of the turbulent vis- 
cosity and the second one is the so-called A-effect due to the 
correlation of the fluctuating velocity components produced 
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Fig. 1 The dynamic turbulent viscosity, computed accord- 
ing to standard mixing-length theory (dotted line) and the 
density (solid line), normalized to their respective values at 
the base of the convection zone, in the model for a solar-like 
star adopted by Lanza (2007). The convection zone extends 
from 0.713 R to the surface, whereas an overshoot layer is 
assumed to extend between 0.673 and 0.713 R with a lin- 
early increasing r) t . 

by the action of the Coriolis force in a rotating star (Riidiger 
1989). 

The conservation of the total angular momentum of the 
convection zone leads to the following boundary conditions: 



6 r = for r = 



R, 



(5) 



where r c is the radius at the lower boundary of the con- 
vection zone and R is the radius of the star. We separate 
the contribution of the turbulent viscosity stress tensor from 
the other terms in Eq. (O and adopt a dynamic turbulent 
viscosity that depends only on the radial co-ordinate, i.e., 
rj t = rjt(r) (see Fig.[T]for details). In such a way, the equa- 
tion for the angular velocity derived from Eqs. (|2), ® an d 
(HJi becomes: 



dn 1 d 
dt pr A dr 
Vt 1 



r rjt 



on 

dr 



d 



pr 2 (1 — p 2 ) dp 



(1-M 2 ) : 



dp 



S. 



where p = cos 6 and the source term S is given by: 

s = - Vr 



pr 2 (\ — p 2 ) 1 
t being a vector whose components are: 

Tj = pr 2 sin 2 #SlM( m )i 

•V; -l{B i B tp + {B' i B' v )) 



(6) 



(7) 



r sin ( 



(8) 



The boundary conditions given by Eq. © can be recast as: 

on 



dr 



for r 



R. 



(9) 



The solution of the angular velocity equation © with the 
boundary conditions (0 can be expressed in terms of an ap- 
propriate Green function (see Lanza 2007 for details): 

n(r,6,t) = ( [V ■T{r',6 , ,t)]G{r,6y,e , )dV' = 
Jv 

= [ (t ■ VG)dV, 
Jv 



(10) 



where the integration with respect to the variables r' and 
9' is extended to the volume of the convection zone V and 
V'G is the gradient of G with respect to the variables r' and 
9'. The second equality in Eq. ( fTOb follows by the appli- 
cation of Gauss's Theorem and the vanishing of r r at the 
boundaries of the convection zone. Note that Eq. ( fTOb is 
valid when the timescale of variation of r is much longer 
than the diffusion timescale of the convection zone t^is — 
(R, — r c ) 2 /v\ , with V\ oc rj t j p the average kinematic turbu- 
lent viscosity of the convection zone. If this assumption is 
not valid, the solution can still be expressed as a linear com- 
bination of terms each depending on an appropriate Green 
function (see Lanza 2007). For the sake of simplicity, how- 
ever, we prefer to work with Eq. ( fTUb . so we retain this as- 
sumption. 

Note that the Green function G depends on r]t and p, but 
does not depend on the angular momentum transport terms. 

We shall obtain information on the averaged intensity of 
the angular momentum transport vector r in the framework 
of two different hypotheses, i.e.: a) the internal angular ve- 
locity is uniform on cylindrical surfaces co-axial with the 
rotation axis (Tay lor- Proudman regime); or, b) the torque is 
localized between, say, r\ < r < r2 within the convection 
zone. 



3 Applications and results 

Let us consider first the case of a stellar convection zone 
rotating in the Taylor-Proudman regime. This is probably 
the case for rapid rotators, like AB Dor and LQ Hya. In 
such a case, Eq. ( fTOb can be written in the form: 

Afi(M) = J Q R C(s')r s (s',t) ds', (11) 

where Ail is the variation of the angular velocity with re- 
spect to a reference rotation state (e.g., that with r = 0), s 
is the distance from the stellar rotation axis in a cylindrical 
polar reference frame and C(s) the lateral area of the cylin- 
drical surface of radius s within the convection zone. If we 
pose: 



M — max 



dG 



ds' 



in the convection zone, 



(12) 



then, Eq. (fTTb can be applied to obtain the inequality (see 
Lanza 2006a for details): 

IAOI 



/ \r s \dV'> 
Jv 



M 



(13) 
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Star 




min {\B S B V \} 




B cq 




AB Dor 


59.6 


0.0385 


0.196 


0.475 


0.23 


AB Dor 


96.7 


0.0624 


0.250 


0.475 


0.60 


LQ Hya 


-48.7 


0.0336 


0.183 


0.495 


0.15 


LQ Hya 


200.9 


0.1386 


0.372 


0.495 


2.52 



Table 1 dil: surface shear in mrad d _1 ; B m [ n : minimum 
magnetic field, obtained for B s = B v ; B cq : mean equipar- 
tition field. Magnetic field intensities are given in Teslas. 



Equation ( TL31 provides us with a lower limit for the aver- 
age radial angular momentum transport term in the convec- 
tion zone. Under the additional hypotheses that the refer- 
ence state corresponds to solid-body rotation with the total 
angular momentum of the star, and r = —(s/£l)B v B, i.e., 
only the azimuthal mean-field Maxwell stresses contribute 
to the deviation of the differential rotation from the refer- 
ence state (cf. Covas, Moss & Tavakol 2005), we find: 

4/3 

mm{\B s B v \} = 8.22x HT 4 (^|) 
/ M \ 2/3 / L \ 1/3 / i? V 5/3 , , , 

and the ratio of the dissipated power to the stellar luminosity 
(Lanza 2006a): 

^.4.40xlO-f^) 4/ VilV /3 
L V 1.75/ \MqJ 

/ t \ -2/3 / p x 4/3 

"(£) (£) " 5 > 

where the surface shear dCl is measured in mrad d 1 , ckml 
is the ratio of the mixing-length to the pressure scale-height, 
M the mass, R the radius and L the luminosity of the star. 
The results of the application of these equations to AB Dor 
and LQ Hya are listed in Table [TJ for «ml = 1-75 and 
adopted stellar parameters M = 1.0 M©, R — 1.05 R©, 
L = 0.618 L© for AB Dor and M = 0.95 M©, R = 0.95 
R© and L = 0.506 L© for LQ Hya, respectively. The first 
column in Table [TJ indicates the name of the star, the sec- 
ond the observed surface shear, as obtained from Zeeman 
Doppler Imaging based on Stokes V data (see Donati et al. 
2003a), the third the minimum Maxwell stress, the fourth 
the minimum mean field obtained for B s = B v , the fifth 
the mean equipartition field derived from the average ki- 
netic energy of convective motions computed according to 
the mixing-length theory, and the last the ratio of the dis- 
sipated power to the stellar luminosity. The minimum field 
strengths turn out to be smaller than the average equipar- 
tition values. However, Zeeman Doppler Imaging provides 
us with an estimate of the radial mean field B s « 0.01 T 
in AB Dor and LQ Hya (Donati et al. 2003b), implying 
that the toroidal mean field is significantly stronger than the 
minimum values listed in Table [TJ i.e., B v w 3 — 14 T. 

www.an-journal.org 



Such strong toroidal fields produce a radial Lorentz force 
that opposes the centrifugal force in a rotating star lead- 
ing to a perturbation of its oblateness and therefore of the 
quadrupole moment of its outer gravitational potential. If 
the active star is a member of a close binary system, this 
leads to a variation of the orbital period through a mecha- 
nism originally proposed by Matese & Whitmire (1983) and 
further discussed by Lanza (2005) and references therein. 

Note that the power dissipated in the turbulent convec- 
tion zone of LQ Hya exceeds the stellar luminosity in the 
case of the largest observed shear, i.e., dYl = 200.9 mrad 
d _1 . This does not pose a serious problem if the duration 
of such episodes of large shear is short in comparison to 
the phases with nearly rigid rotation because the large en- 
ergy reservoir represented by the thermal content of the stel- 
lar convection zone will average out the energy loss over a 
timescale comparable with the Kelvin-Helmoltz timescale 
of the envelope. Moreover, the uncertainty in the determi- 
nation of the largest dVl excursion of LQ Hya (cf. Donati et 
al. 2003a) makes the estimate of the corresponding energy 
dissipation rate largely uncertain. 

Let us now consider the case in which magnetic stresses 
are localized close to the base of the convection zone, with- 
out assuming a specific internal rotation regime. In this case, 
an expression similar to Eq. ( TT31 can be obtained (cf. Lanza 
2007). Specifically, we assume that r = — (r sm9/p)BB v 
is localized in the overshoot layer, i.e., 
0.67 < (r/R) < 0.71 and that the differential rotation ob- 
served at the surface is equal to that at r = 0.99i?. 

In the case of the largest shear observed in LQ Hya, this 
yields a minimum average Maxwell stress over the over- 
shoot layer {B r B v } . ~ 0.063 T 2 , which gives fields of 
the same order of magnitude of those obtained above under 
the Taylor-Proudman hypothesis. 

4 Discussion and conclusions 

We modelled the observed time variation of the differen- 
tial rotation in AB Dor and LQ Hya under the hypothe- 
ses that the azimuthal Maxwell stresses rule the changes 
of their surface shear and the internal angular velocity de- 
pends only on the distance from the rotation axis (Taylor- 
Proudman regime). We obtained that the average intensity 
of the mean field Maxwell stresses is \B S B V \ ~ 0.03 — 0.14 
T 2 , implying azimuthal mean fields B v ~ (3 — 10) T for 
-B s ~ 0.01 T. Similar Maxwell stresses are obtained if the 
magnetic torque is assumed to be confined within the over- 
shoot layer 0.67 < (r/R) < 0.71 and no restrictions are 
imposed on the internal rotation law. 

It is interesting to note that azimuthal magnetic fields of 
3 — 10 T, occupying a sizeable fraction of the convection 
zone, have been invoked to explain orbital period changes 
observed in late-type active binaries (Lanza, Rodono 
& Rosner 1998; Lanza & Rodono 2004; Lanza 2005, 
2006b). An a-effect related to an instability of the mag- 
netic field itself (e.g., magnetic buoyancy instability, Bran- 
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denburg & Schmitt 1998; or magneto-rotational instabil- 
ity, Riidiger et al. 2007) seems to be necessary to produce 
such super-equipartition fields, possibly acting in combina- 
tion with differential rotation. The energy dissipated by tur- 
bulence, estimated according to standard mixing-length ar- 
guments, may exceed stellar luminosity in the case of the 
largest surface shear observed in LQ Hya. However, the 
thermal equilibrium of the convection zone can be signif- 
icantly affected only if those large shear episodes last more 
than ~ 10 — 20 per cent of the time. Note also that a mixing- 
length estimate for the turbulent viscosity may not be ap- 
propriate for a rapidly rotating star (see Lanza 2006a for 
details). 

In the present work, we adopted a point of view analo- 
gous to that of Covas et al. (2005) who modelled the varia- 
tion of stellar differential rotation considering only the 
torque exerted by the Lorentz force and neglecting the roles 
of meridional flow and A-effect. As a matter of fact, it is 
difficult to evaluate the perturbations of the meridional flow 
and of the Reynolds stresses produced by the magnetic field 
because they depend critically on the approximations made 
in the treatment of stellar turbulence in a rotating star. Nev- 
ertheless, alternative models for the variation of differential 
rotation have been investigated, such as those based on a 
time-dependent component of the meridional flow (Rem- 
pel 2006, 2007) or the magnetic quenching of the A-effect 
(Riidiger et al. 1986). The present approach can be gener- 
alized to obtain amplitudes of the perturbations of the cor- 
responding terms in Eq. ([8]), but we shall not pursue this 
application here (see Lanza 2007). 

Finally, it is interesting to note that some inference on 
the amplitude of variation of the surface shear in the case 
of very active stars can be obtained not only by means of 
Doppler imaging techniques, but also by an appropriate ana- 
lysis of their long-term wide-band photometry (e.g., 
Rodono et al. 2001; Messina & Guinan 2003). For example, 
in the case of LQ Hya, it is worth comparing the Doppler 
imaging results by Donati et al. (2003a) with the photomet- 
rically determined seasonal rotation periods by Kovari et al. 
(2004). 
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